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To analytically derive the field of circular polarized incident laser, it is nec- 
essary to note that in the rotated coordinations, Eq. (1) still holds with phase 
mask function e''=*('*'''^') = 1 for the excitation beam and e*'''*^^'''^') = e'^'^'+f ) 
for the depletion beam. It yields 

eL'(-P) = -iAiy{Io - /2COs2(/)p) 
eiy'(-P) = iMyh sin2(/)p 
eL'(-P) = -2Aiyhsm(j)p 



for excitation and 



A2y[yi cxp{i(j)p) + ij^2cxp(-i0p) - i^3exp(i30p)] 

— ^ [J^2 exp{-i4>p) + J^3 exp(i30p)] 
-A2y[^4 + Acxp(i20p)] 



for depletion. 
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In Eq. (4) the parameters are defined as 
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= &41 + bi2{k2Pp?' 
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= h{k2Ppf' 



(SI) 



where oqi = [B(s, 1,3/4) + B(s, 1, 5/4)]/2,ao2 = [S(s, 2, 3/4) + B(s, 2, 5/4)]/8, 
ai = S(s,2,3/4)/4, 03 = [i?(s, 2, |) - i?(s, 2, |)]/16, W = [i?(s, 3/2, 3/4) + 
B(s,3/2,5/4)]/4, 62 = [5(5, 3/2, 3/4) -B(s, 3/2,5/4)1/4, 641 = B(s, 3/2, 3/4), 
642 = B{s, 5/2, 3/4)/4, 65 = B{s, 5/2, 3/4)/8. 

To obtain Eq.(5) it is important to keep in mind that without phase plate 
Idep = D^clqi- If we defined the depletion power as when peak fluoresce 
intensity into half, we also has Is = D^a^^. So depletion later power D^/D^ 
in Eq.(4) is equivalent to peak intensity ratio Idep/Is in Eq.(5). For continuous 

wave, Iex{p)v{Idepip)) = Iexi0)/2, r]{Idep) = 1/(1 + Idep/Is): lex = C^[ali + 

2(aoiao2 + a^)(fcip)2]and Idepip) = Dlbl{k2p)^ yield 

FWHM = 

V-4(aoiaQ2+af) L bj / A,^ p /step ^^'^^ 



aoi V 2(aoiao2+al) ^ ' 

And for continuous wave, Iex{p)r]{Idepip)) = Iexi0)/2, 7?(/rfep) = exp{-ln2Idep/Is), 
lex = C^al^expl + 2(aoiao2 + a?)/aoi('=iP)^ and /rfep(p) = Dfbl{k2p)^ yield 

FWHM = .QON 

/ -2(Q0ia02+Qj) /-, I (i»2)b; . Aex \2 JsTiSD V'^'Jj 

Then in Eq.(5) factor a and (3 can be obtained. 

To verify our code STED3D, first we simulated the excitation intensity on 
xy and xz plane in Fig. SI. The FWHM is 256.8 nm and 605.7 nm respectively, 
in good agreement with result from PSFlab [19], thus verifying our code. Then 
we plotted in Fig. S2 the intensity on the xy plane of the focus for the depletion 
beam with left-circular polarized (a) , linear polarized (b) , and right-circular po- 
larized laser (c). It confirms our conclusion that only the left-circular polarized 
depletion beam results in zero intensity in the focus. Also we plotted in Fig. 
S3 intensity of left-circular polarized depletion beam with Eq. (S2) in (b). The 
simulated STED distribution in both linear and circular polarization agreed well 
with the results of [23], and again the matching results verified our derivation. 

In the discussion part we referred the saturation power corresponds 
to about twice the saturation intensity Is- To illustrate this, we plotted the 
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Figure SI: Excitation intensity on (a) xy plane and (b) xz plane. 
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Figure S2: Intensity on the xy plane of the focus for the depletion laser with 
(a)left-circular polarized, (b) linear polarized, and (c) right-circular polarized 
laser. 



nonlinear curves of integration depleted fluorescence efficiency and depletion 
intensity in Figure S4. And we also listed different factor ^ with various depletion 
wavelengths of STED dyes in Table SI. 
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Figure S3: Intensity of the left-circular polarized depletion beam with Eq. (1) 
in (a) and Eq. (2) in (b) 
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Figure S4: Integration depleted fluorescence efhciency and depletion intensity. 
Simulation results are plot in blue solid for CW-STED and red dash for pulsed- 
STED. Then fractal (blown dot) and exponential (green dash and dot) fits are 
performed for two conditions respectively. 



4 



Dye 


Wavelengths /nm 


Pulsed STED 


CW STED 


Excitation 


Depletion 


X 


X 


ATTO 532 
Alexa Fluor 488, Chrorrreo 488 
FITC, DY-495, Citrine 
GPP 
Citrine, YPP 
Pyridine 2, RH414 
Alexa 594, Dylight 594, ATTO 594 
ATTO 633, ATTO 647N, 
Abbcrior STAR635 


470 

488 

490 
490 
554 
570 

635 


603 

592 

575 
600 
745 
690 

750 


0.31 

0.36 

0.42 
0.35 
0.26 
0.36 

0.42 


0.40 

0.45 

0.53 
0.44 
0.33 
0.45 

0.52 


Table SI: Typical values 


of (3 for 


different 


STED sys- 





terns. Dates of excitation and depletion wavelengths are from 
http : //nanobiophotonics . mpibpc . mpg . de/ old/dyes/. Numerical aper- 
ture of the lens is set to be NA = 1.4, and a = 2.15 for pulsed STED, a = 1.82 
for CW STED. 
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Abstract 

Stimulated emission depletion (STED) can achieve optical super-resolution, 
with the optical diffraction limit broken by the suppression on the periphery 
of the fluorescent focal spot. Previously, it is generally experimentally accepted 
that there exists an inverse square root relationship with the STED power and 
the resolution, yet without strict analytical description. In this paper, we have 
analytically verified the relationship between the STED power and the achiev- 
able resolution from vector optical theory for the widely used 0-27r vortex phase 
modulation. Electromagnetic fields of the focal region of a high numerical aper- 
ture objective are calculated and approximated into polynomials, and analytical 
expression of resolution as a function of the STED intensity has been derived. 
As a result, the resolution can be estimated directly from the measurement of 
the saturation power of the dye and the STED power applied. 
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Introduction 



In 1873, Ernst Abbe derived the diffraction limit of conventional microscopy 
as dAbbe = X/2NA, which is also regarded as the resolution limit for optical 
microscopy [1]. It was not until the past few decades with the development of 
super- resolution microscopy techniques that this limit has been surpassed [2]. 
One of these techniques, stimulated emission depletion (STED), belongs to the 
optical nanoscopy family of reversible saturable optical fluorescence transitions 
(RESOLFT) [3], in which the point spread function (PSF) of the fluorescence is 
engineered to obtain a higher resolution which breaks the optical diffraction limit 
[4]. In particular, a spatially modulated depletion beam is spatially overlapped 
on the excitation focal spot to effectively remove the peripheral fluorescence 
distribution with either STED [5] , or by exciting the electron to a triplet state 
(ground state depletion, GSD) [6], or through excitation state absorption (ESA) 
to excite the electron to a higher energy level [7] . To achieve super-resolution in 
the lateral plane, a widely used method is to place a 0-27r vortex phase modulator 
in the optical path of the circular polarized depiction beam [8, 9, 10, 11]. 

To estimate the theoretical super-resolution power of STED, the most widely 
used expression is the square root law in which the full-width at half-maximum 
FWHM = ho/y/l + I/Is which was first developed by Westphal and Hell [12] 
and later by Harke et al. [13]. Although partially explains the experimen- 
tal results, it is an empirical theory without solid theoretical foundation, as 
the depletion intensity distributions discussed in previous work are not derived 
from the diffraction theory. The scalar solution to the depletion intensity of a 
0-27r phase plate has been derived by Watanbe et al. [14]. Yet without con- 
sidering the vectorial nature of light, such theory can only be applied with a 
small aperture objective where the paraxial approximation is applicable. As a 
result, this theory fails for objectives with large apertures. On the other hand, 
it is possible to numerically simulate the intensity distribution at the focal spot, 
without analytical expression [15, 23]. So far, despite the huge progress in the 
experimental instrumentation and application of STED in optical nanoscopy, 
accurate estimation of the performance of a STED system under a known de- 
pletion power is still lacking. It is therefore of great importance to establish 
a new formula describing the maximum achievable resolution of STED optical 
nanoscopy, based on the relationship between the power of the depiction beam 
and the measured power when the detected fluorescence intensity is depicted to 
its half. 

In this paper we use a vectorial theory to investigate the laser field in the 
neighborhood of focus. We present the diffraction analysis in Results, where 
analytical excitation and depletion light intensity distributions are given and 
then approximated with polynomials. It is demonstrated that only the chirality 
of the circular-polarization matching the direction of the 0-2tt vortex phase plate 
can generate a central zero doughnut intensity distribution. Furthermore, with 
a second-order approximation, expressions of FWHM with different saturation 
functions are derived from system parameters in both CW-STED and pulsed- 
STED. The previous resolution law FWHM = h^l \/\ + ///, is validated, and 
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the relationship between the depletion power and the half fluorescence inhibition 
power are developed. Then we numerically simulate the electromagnetic fields, 
and compare its results with our analytical results. Finally wc discuss the error 
induced by our approximation, and the scope and limitation of the application 
to this theory. 



Results 

Wolf et al. derived the vectorial theory of diffraction under asymptotic and 
the optical rays approximation [17, 18]. For a linear polarized excitation laser 
beam, the electric beam in objective space can be written as: 



TT 



cos5 6lsme'{cos 6* + (1 + COS 61) sin^ 0}e''''''' ™" ' 



Jo 

a ^27r 



e=— / / cos^6lsin6'(l-cos6l)sin0cos0e''=''''™"'d6lfi0 



iA 

TT 



^0 
a p27r 



Jo 



cos^ esin^ecosche''"''' """^ 



} (1) 



The coordinates arc illustrated in Fig. ??. In an idealistic Gaussian image 
system, the above expression yields: 

exiP) = -iAix{Io + l2Cos2cj)p) 
ey{P) = —iAixl2sin24)p 
eziP) = -2Ai.Jicoscj)p 

where Aix = kifeix/2, ki = 2'K/Xex and eix is the linear excitation electric 
field along the x-axis, and on focal plane integration /q, /i, I2 are given by 

/•a 

10 ~ cos"^ 8 smd{l + cos9)Jo{kippSm9) eyii>{ikiZpCOs9)d9 

Jo 

/•a 

11 = cos^ 9 sin^ 9Ji(kippSm9) cxp{ikiZpCos9)d9 

Jo 

/•a 

1-2— cos^ 9 sm9{l — cos9)J2{kippSm9) exp{ikiZpCos9)d9 
Jo 

where pp = {Xp + UpY^^- In STED a 0-2tt phase mask is generally used, as it 
can generate the most efficient inhibition pattern for fluorescence depletion [3] . 
The electric field of a depletion beam with distinct polarization states can be 
calculated [15]. Notice /g^'^ expjmc/) + p cos 0}(i(/> = 27ri"J„(p), and substitute 
exp(ifcrpcose) = e*''"^J'™"^e''=''f ^^^^('^-'^p^e*''*^^''^) into Eq. (1), where the 
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phase function of 0-27r phase plate is e 



.ik<s>{0.,,l>) ^ gi0^ derive: 



y 



Mx\-^\ exp{i(j)p) - - J^2 cxp{~i(j)p) + -J^3 cxp(i3(^p)] 
^ [=^2 exp(-z0p) + J^3 exp{i3(j)p)] 



> 



(2) 



= J^2x[^4 - -<^5 exp(i2(?!)p)] 

Here = fe/e2a;/2 is proportional to linear depletion electric field, and fc2 = 
27r/A(ie- The angular integrations ^1 — J^s are defined as 



Comparison of the 3D intensity profile with result of PSF Lab [f 9] can be 
seen in Fig. SI. 

For circular polarization incidence, the electric field can be obtained by co- 
ordinate rotation [20]. The electromagnetic field now consists of contributions 
of two linear polarized incident fields along the x and y axis: 



where = Tp;(^'p ~ Sp;4>'p ~ 4>p ~ -I.Thcn the intensity of laser field could 
be calculated by / = E^E* + EyE* + E^E* for both the excitation and the 
depletion beam, and residual intensity is Iresiduai — lex ^vi^de), where ri{Ide) = 
1/(1 + Ide/Is) for continuous beams and rj{Ide) = exp[— (In 2)/(ie//s] for pulsed 
STED [21]. 

Now consider the light intensity distribution on focal plane Zp = 0. Since 
Jy(x) = k\T(k+u+i) (f )^^~'"''' the approximation to the second-order can be 
used: Jq{x) = 1 — Ji{x) ~ |, and Ji{x) = Define s — sin^ a and 



B{s,p,q) = JotP-'^il - ty-^dt = 2£sm^P'^6cos^i-^6d6, which is the in- 



complete Beta function, we can express Ii — ^5 as polynomials of pp with their 
coefficients given by aoi — &5, and their coefficients are incomplete Beta functions 
of s, as shown in Supplementary Information. For a linear or elliptically polar- 
ized laser, it is always possible to rotate coordinates to make the polarization 




Ex{rp,ep,(j)p) 
Ey{rp,Op,(j)p) 
Ez{rp, Bp, 4>p) 



exirp,9p,(l)p) - e'y,{r'p,e'p,(t)'py 
ey{rp,ep,(l)p) + e'.^,{r'^,e'p,(l)'p) > 
ez{rp,0p,(l)p) + e^,{r 0' (f)' ) , 



(3) 
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a standard ellipse with Ai^ = Ci, Aiy = iC2, A^x = Di, A2y = iD2, where 
Ci,C2,Di,D2 are real. Then the excitation and STED doughnut intensities can 
be expressed as: 

leApp, 0p) = {Cf + Cl)al^ + 2[{Cl + C|)(aoiao2 + a?)' 

IdepiPp, (t^p) = {Dl - D2fbl, + {Dl + Dl)bl{k2Pj>f I (4) 

+ [Dl - D2f{hjl + 2biihi2){k2Ppf 

- (6162 + 264165) (I?? - Dl)cos2(l)p{k2Ppf , 

This expression gives an analytical description of the excitation and depiction 
intensities in the vicinity of the focus. A special case is that the depletion beam is 
left-hand circularly polarized (the same as the direction of the vortex direction) 
with Ci = C2 = -^C and Di = D2 = -^D, which indicates that the center 
of donut has zero intensity. Then I^x = C^i^^oi + 2(aoiao2 + al){kip)^] and 
Idep = D\h\{k2pY ■ For an elliptically polarized laser or linearly polarized laser 
[22], Dl 7^ D2 leads to Ide 7^ at the focus. This is derived from second-order 
approximation, but still holds true for a strict analytical solution because all 
the higher-orders annihilate at ?' = 0. The results obtained with our methods 
can also be validated with [23], as shown in Fig. S2. The full-width at half- 
maximum (FWHM) of STED PSF could be estimated accordingly, with the 
decay functions for CW or pulse cases described in SI. 



Discussion 

The intensity Ig could be measured as the following procedures: (1) measure the 
peak fluorescence intensity of excitation beam (here the size of the fluorophore 
is assumed to be very small, i. e. a delta function); (2) introduce the depiction 
beam without phase modulation, whose PSF should be overlapped with the 
PSF of the excitation; (3) adjust the depletion power until the peak fluorescence 
intensity is depleted to its half [13]. As a result, the peak depletion intensity 
is defined as Ig- When the CW saturation functions rj{Ide) = 1/(1 + Ide/ Is) is 
considered, to its second-order [21]: 

FWHM = , ^^"^ f5) 

aNA^l + Phep/Is ^ ' 

where Idep is the non-modulated peak intensity of the depletion beam focus, 

and a = TT^-A{a^ia^2 + a?)/aois, /? = 2(ao.al+aj) (l^)'- For NA = 1.4 
and n = 1.5, Xex = 635nm, Xde = 760nm, this expression yields FWHM = 
Xex/{2.15NA^yi + 0.168Idep/Is)- For pulsed-STED with saturation function 
77(1) = exp(— Zn2///s), it is more convenient to rewrite excitation into an ex- 
ponential function [21] and thus a = tt-^/— 2(aoiao2 + a^)/(/n2)/aois, /3 = 
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■5y(j^)^. With the same system parameters, this expression yields 



FWHM = Aex(l + 0.233/dep//s)"^/Vl-83AfA. We found the inverse root law 
still holds in vectorial theory except for modification factors a and /3, which are 
incomplete beta functions with aperture angle as their parameters. The rela- 
tionship between the depletion intensity ratio I / Is and the fluorescence intensity 
ratio Pdet/Po is shown in Fig. S4. 

In experiment, the measurement of intensity is often replaced with measuring 
of the laser power, and then calculating the intensity distribution for modulated 
and unmodulated STED PSFs [24]. Stimulated emission intensity Is is often 
defined as the intensity of fluorescence dropped into its half [13, 21], or 1/e 
[12, 25]. The resolution of STED is often described with the relationship of 
PjSPsi where P is the depletion power and !^s is defined by the depletion 
power when the detected fluorescence intensity is reduced by half, if we overlap 
the non-phase modulated depiction PSF with the excitation PSF [24, 13]. 

Wc used STED3D to simulate the electromagnetic field distributions around 
the focus. In Fig. we plotted the excitation (a) and the depletion (b) intensities 
on X-axis of the focal plane. In (a) the simulation intensity was plotted in red, the 
polynomial approximation of Eq. (4) was plotted in blue, and the exponential 
approximation was plotted in green. In Fig. (b) the numerical solution to the 
rigid EM field intensity was plotted in red, and polynomial approximation of 
Eq. (4) was plotted in blue. As expected, approximations agreed well with the 
numerical solution of the EM field in the vicinity of focus, but diverged with the 
increase of the radius. Then we plotted in Fig. the simulated FWHM (red) 
and the calculated results in Eq. (5) (blue). Figure (a) is the CW case with 
fraction saturation function and Figure (b) is the pulsed case with exponential 
saturation function [21]. In both figures the two curves overlap better with the 
increase of depletion intensity. 

In above discussion we followed Wolf's diffraction theory and Harke's ap- 
proximation to derive an analytical resolution expression. Yet, it is difficult to 
measure the focal intensity of a high-NA objective. Practically, the intensity is 
calculated with the measurement of the laser power and diffraction distribution 
[13, 24, 26]. First the total fluorescence power with excitation beam is record by 
a confocal point detector [13, 24], and then depletion beam is introduced with- 
out phase modulation. Reduced detected fluorescence power is recorded with 
respect of the depiction power, and the exponential-fitted saturation power is 
defined as S^s [26]. Therefore, it is of great importance to develop relationship 
between the STED resolution and the STED power [26]. Numerical evaluation 
of the factor reveals P — i^s approximately corresponds to / = lis (Appendix 
Fig. 6). Then the resolution vs. STED power relationship becomes [24, 27] 



where P is the depletion power, ^s is the saturation power when the fiuo- 
rescence is reduced by half. Our representation of resolution then becomes 



FWHM = Ae^/(2.15iVyl^l + QAP/^s) and FWHM = \ex / {l.^2N A^l + Q.bP/^s) 



FWHM = 



aNA^l + xP/^s 



(6) 
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for CW-STED and pulsed-STED, respectively. As an example typical values of 
parameter x for different dyes are listed in Supplementary Information Tab. SI. 

In this paper we have applied a polynomial approximation of the Bcsscl in- 
tegrals, to yield an analytical expression of STED with vectorial waves. The re- 
sults were validated with previous reported numerical simulations of the electro- 
magnetic wave. The relationship of 1/ y^l -I- I /Is has been validated in both CW 
and pulsed STED cases. We have given new sets of resolution estimation equa- 
tions, to predict the resolution achievable for a STED system: the expression 
of resolution, incident depletion intensity and half-peak-fluorescence depletion 
intensity is analytically derived with incomplete Beta functions, and the rela- 
tionship of half-peak-fluorescence depletion intensity and half-total-fluorescence 
depletion power was established by simulation. The result can be extended to 
other RESOLFT type super-resolution techniques where the fluorescence point 
spread function is modulated with another doughnut shape point spread func- 
tion through 0-27r phase modulation. 

In the derivation of Eq. (5), asymptotic expansion at ?' = makes the 
approximation only available in the vicinity of the focus. In Fig. we found 
that the curves of analytical and simulation laser intensity overlap well when r 
was small but diverged with the increase of r. In Fig. we found that when 
r < SOnm, the approximated excitation and depletion intensity will induce 
errors of 5% and 8%, respectively. As a result, our theory applies to high 
aperture objective lenses with NA > 1, or those with high depletion intensity, 
roughly when the FWHM is less than A/3. For smaller aperture objective lens, 
Watanabe's theory [14] with scalar diffraction theory [28, 29] can provide an 
acceptable estimation. For the same reason, in Fig. the analytical FWHM 
curves better predicts results with error less than 5% when / > lO/s. With 
the development of the high efficiency fluorescent dye and through application 
of larger STED power. STED resolution of 6 nm has been demonstrated [4], 
and typical STED system resolution is < 80nm. As usually the best STED 
nanoscopy resolution at a certain depletion power is concerned, rather than the 
resolution at moderate power, our result has a wide scope of application in the 
accurate estimation of the STED (or RESOLFT) system performance. 

There are several improvements that can be made in the future. One is 
introducing the Gaussian spatial distribution of the incident lasers, which can 
be done by inserting a Gaussian factor in the integration of Iq — J^s. Since the 
beam is generally expanded for overfilling the back aperture of the objective, the 
planar wave can serve as a very good approximation to the experimental situa- 
tion. Another possible improvement would result from retaining higher orders 
in the Bessel summation. It may reduce the error when the radius increases, at 
the cost of much complexity in expression. 
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Figure captions 



Fig. 1 Schematic diagram of the diffraction of a high NA microscopic objective. 
Fig. 2 The intensity profile of (a) excitation and (b) depletion beams on x-axis 
of the focal plane. Simulation intensities are plotted in red solid, polynomial 
approximation of Eq. (4) arc plotted in blue dash, and the exponential approx- 
imation is plot in green in (a). 

Fig. 3 The relationship between the depletion intensity and its resulted res- 
olution. Simulation results are plotted in red solid for CW-STED in (a) and 
pulsed-STED (b). Then predicted resolutions of Eq. (5) are plotted in blue 
dash in two figures, respectively. 
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